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We propose a physical mechanism which enables permanent Rabi oscillations in driven-dissipative 
condensates of exciton-polaritons in semiconductor microcavities subjected to external magnetic 
fields. The method is based on incoherent excitonic reservoir engineering. We demonstrate that 
permanent non-decaying oscillations may appear due to the parity-time (PT) symmetry of the 
coupled exciton-photon system realised in a specific regime of pumping to the exciton state and 
depletion of the reservoir. For effective non-zero exciton-photon detuning, permanent Rabi oscilla¬ 
tions occur with unequal amplitudes of exciton and photon components. Our predictions pave way 
to realisation of integrated circuits based on exciton-polariton condensates. 


Introduction. —Nowadays, semiconductor microcavi¬ 
ties with embedded quantum wells serve as a solid-state 
laboratory for fundamental studies of dynamic, coher¬ 
ent, nonlinear and quantum effects in nonequilibrium 
ensembles of bosonic quasiparticles: exciton-polaritons 
mm- Various approaches to photonic information pro¬ 
cessing with the use of microcavity polaritons have been 
proposed, see, e.g., [3H5]. In particular, polariton cir¬ 
cuits, or neurons, are proposed in [6j and then experi¬ 
mentally examined by Ballarini et al. in [7] at extremely 
low pump intensities. Although nonequilibrium exciton- 
polariton Bose-Einstein condensation has been already 
observed in many labs, the requirement to use high Q- 
factor microcavities with high reflectivity Bragg mirrors 
limits wide application of such phenomena for practical 
purposes. Even in highest quality microcavities used in 
current experiments the leakage of photons leads to dis¬ 
sipation effects occurring on the scale of tens of picosec¬ 
onds. In particular, this leads to a rapid attenuation 
of Rabi oscillations occurring between exciton and pho¬ 
ton components of a polariton state, which represent a 
particular interest for quantum information applications. 
Recently, some of us have proposed [5] a way to improve 
the coherence time of polariton Rabi oscillations by stim¬ 
ulated pumping of a Rabi-oscillator from a permanent 
thermal reservoir of polaritons. The increase of the the 
coherence time of polariton Rabi-oscillations in the pres¬ 
ence of cw pumping has been experimentally observed [5], 
and manifestations of this effect in spatial dynamics of 
exciton-polaritons have been revealed experimentally [2] 
and analyzed theoretically [101. Recently, Voronova et al 
m has studied theoretically the non-linear regime of po¬ 
lariton Rabi-oscillations and their interplay with exciton- 
photon Josephson oscillations. 

In the Letter, we demonstrate theoretically realisa¬ 
tion of permanent Rabi oscillations between excitonic 
and photonic components of a driven-dissipative spinor 
exciton-polariton condensate in semiconductor microcav¬ 


ities at specific pumping conditions. We study the regime 
where the microcavity system obeys parity-time (PT) 
symmetry conditions. Originally, the PT-symmetry ap¬ 
proach has been proposed by Bender et al pj} [13] to 
demonstrate that non-Hermitian Hamiltonians can pos¬ 
sess entirely real energy eigenvalue spectra. In the past 
decade, PT-symmetry features have been demonstrated 
for various systems in photonics [14H19I , condensed mat¬ 
ter physics [20] and metamaterials [21]. Moreover, PT- 
symmetry approach has been extended to electronic cir¬ 
cuits |22] which makes it highly relevant to polariton- 
ics. Physically, PT-symmetry requirements can be easily 
understood by using the system of two linearly coupled 
oscillators (dimers, or waveguides in photonics [171 fT3]). 
The system obeys PT-symmetry criteria if the rate of dis¬ 
sipation for one of the oscillators is exactly equal to the 
rate of gain in the other oscillator. This condition can 
be achieved in a system of two coupled polariton conden¬ 
sates with distributed dissipation rates [23], that results 
in the conservative Hamiltonian dynamics. It is impor¬ 
tant to note that even if the exact compensation between 
gain and losses in a two-nrode system is not achieved, the 
concept of quasi-PT-symmetry can be introduced 231 . 
Here we demonstrate that the regime of permanent Rabi 
oscillations can be realized for the exciton-photon sys¬ 
tem in a microcavity, and the condition of PT-symmetry 
may be achieved if the gain in the excitonic component 
is compensated by losses from the photonic component. 

Basic equations .—We consider a coupled exciton- 
photon system in the presence of the external magnetic 
field that produces the Zeeman splitting hAz of exci¬ 
ton levels, and an incoherent excitonic reservoir, pumped 
by cw-field P. This non-resonant pumping can be real¬ 
ized both optically or by electronic current injection [251 . 
We assume that the wave function of driven-dissipative 
exciton-polariton condensate consists of photonic and ex¬ 
citonic components, (f>± and %±- We will neglect the pos¬ 
sible spatial degrees of freedom in what follows, assuming 
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the condensate to be at zero momentum state. In this 
case, both components of the condensate and the reser¬ 
voir can be described by Boltzmann kinetic equations: 


X± = 2 (Px[N±\ - 7x) X± + i$±X± 

<t>± = - in x±> 

N±=P- 7 flIV± -Px[AI±] |x±1 2 , 


(l a ) 

(lb) 

(lc) 


where dots denote time derivatives and 2f2 is a polariton 
Rabi splitting frequency. The subscript “+” (“—”) cor¬ 
responds to the spin projection parallel (antiparallel) to 
the vector of magnetic field. In Eqs. Q <5± = A ± A z 
are effective photon-exciton detunings determined by the 
detuning A = up — u>x of the cavity mode and exciton 
frequencies in the absence of magnetic field and by the 
Zeeman splitting A z\ lx, 7 p and 7# are the exciton, cav¬ 
ity and reservoir damping rates, respectively. In Eqs. ([T]) 
the term containing px [IV±] describes the pumping of the 
exciton state by stimulated scattering from the reservoir. 
We shall consider two possible scattering mechanisms de¬ 
scribed by pumping terms p\ and p 2 '- 

Px [1V±] = pi = RiN ±, (2a) 

px[N ± }=p 2 = R 2 Nl\ X ±\ 2 . (2b) 


The Eq. (2a) implies the acoustic phonon assisted 
pumping with the rate R\ , see [26], while Eq. (2b) de¬ 
scribes the exciton-exciton scattering where excitons pos¬ 
sessing momenta — k and k scatter into the condensate 
state with the momentum k = 0. Note that in both 
cases, the scattering feeds both upper and lower exciton- 
polariton branches. Our model accounts for the decay of 
both exciton and photon components of the condensate 
but neglects incoherent processes that lead to relaxation 
between upper and lower polariton branches. 

To start with, let us neglect scattering of excitons with 
opposite spins so that Eqs. 0 can be solved for each of 
spin components separately. We first examine the sub¬ 
system with spin component omitting subscript for 
simplicity. 

The dynamics of the exciton-polariton system demon¬ 
strates a threshold behavior at the specific value P th of 
the cw-pump which is dependent on the particular mech¬ 
anism of the exciton pumping px [IV] [26] .At P < P th 
the exciton pumping rate px[N] is comparable to or 
much smaller than dumping rates 'yx and 7p which, in 
turn, are much smaller than the Rabi splitting Q (see, 
e.g., [27]). In this case one can assume the reservoir pop¬ 
ulation N and the pumping rate px[N] to be time in¬ 
dependent on a time scale of the Rabi oscillation period 
~ I2 -1 . In particular, for the pumping term given by 
Eq. (2a I, eliminating Eq. 0 ) we find eigenfrequencies of 


the exciton-photon system 

i . S 

U 1,2 = 4 (7 -Px) + 2 

± \ \j 4n2+ i 5 ^ ^2 ^ P ~ 1X +Px ^) ’ 

characterizing steady-state solutions x = Xo e “* and 
(j> = 4>o e lult - Physically, Eqs. ([3]) determine frequencies 
of upper (wi) and lower (uj 2 ) polariton branches, cf. (27j. 
In Eq. ([3]) we denoted 7 = 7 p + 'yx- In this case, from 
Eq. 0) one obtains P{ h = j^r/Ri- 

The PT-symmetry of the system 0 appears if the her- 
mitian part of Hamiltonian is symmetric and the anti- 
hermitian part is anti-symmetric with respect to permu¬ 
tation of excitonic and photonic components. The first 
condition leads to 


(5 = 0, (4a) 

and the second condition requires that the losses from 
the photonic component should be fully compensated by 
the gain of the exciton one, 


Px - lx = 7 p- 


(4b) 


Clearly, if these conditions are satisfied the eigenfrequen¬ 
cies 0 become fully real. We note that in this case 
the amplitudes of oscillations of exciton and photon pop¬ 
ulations become equal |y| 2 = \(f>\ 2 = \A\ 2 . Conditions 
(0,b) can be experimentally realized by tuning of the 
pumping of the excitonic component with variation of 
the cw pumping strength. 

Figure [l] reveals the decay dynamics of a coupled 
matter-light system determined by Eqs. (0 for different 
pumping mechanisms. The chosen parameters are rel¬ 
evant to experimentally accessible GaAs semiconductor 
microstructures, cf. [T(. In the Letter for simplicity we 
represented excitonic |x| 2 and photonic | <f>\ 2 populations 
in the units of initial condensate density. As it is seen 
from Fig. [I] and follows from Eq. (0, the characteristic 
time tr a 2/ (7 — px) of Rabi oscillations decay can be 
increased by manipulating the pump field parameter px ■ 

Condition 0) is a criterion for realisation of perma¬ 
nent Rabi oscillations {tr —> 00). It is easy to show 
that such a regime cannot be achieved below threshold, 
at P < P th , see Fig. [I] Actually, since pumping rate 
px = px[N(t)] depends on the number of particles N(t) 
in the reservoir, it varies as a function of time, and the 
condition (0) just cannot be preserved in the whole time 
interval. This indicates that the features of the combined 
exciton-photon system are strongly affected by the reser¬ 
voir dynamics, as the inset in Fig. [t] shows. 

In the case of exciton pumping predominantly gov¬ 


erned by the exciton-exciton scattering process Eq. (2b), 


the dynamics of the system is more complicated. The 
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FIG. 1. (Color online) Time dependence of the exciton pop¬ 
ulation \x\ 2 for 6 = 0 and for exciton state pumping pi 
(magenta curve) and p 2 (yellow curve) for cw-pump P be¬ 
low threshold P th . Dash-dotted black curve indicates decay 
dynamics in the absence of reservoir. On the inset the rele¬ 
vant population of the reservoir N for exciton state pumping 
Pi (magenta dashed) and P 2 (orange solid) is shown. Hor¬ 
izontal lines correspond to the stationary reservoir popula¬ 
tion, i.e. N = P/'iR. Parameters are: ■yx = 0.01 ps -1 , 
7 p = 0.1 ps” 1 , htl = 2.5 meV, 7 _r = 0.003 ps” 1 , hRi = 
0.01 pm 2 ps _1 and hR 2 = 0.001 /rm 6 ps _1 . Pumping rates are 
Pi — 0.02 /rm” 2 ps” 1 and P 2 = 0.035 /rm” 2 ps _1 - for ma¬ 
genta and yellow curves, respectively. Initial conditions are: 
x(0) = 0, m = 1, N( 0) = P/ 7 fl. 


threshold value P.j h can not be derived directly from 
Eq. (§ because of the dependence of P 2 term on the pop¬ 
ulation |x| 2 of the excitonic component. Figure [l] (yellow 
curve) demonstrates the behaviour of exciton and pho¬ 
ton components below threshold P < P 2 1 , that was de¬ 
termined by numerical simulations. Although both com¬ 
ponents of the condensate are initially amplified, then 
the population of |\| 2 (and hence the pumping term 
P 2 = R 2 N 2 |xl 2 ) attenuate very rapidly. 

Permanent oscillations. —Now we analyze permanent 
Rabi oscillations regime, that occurs above threshold 
p pth an( ^ can kg p reservec i as long as the PT- 
symmetry condition is fulfilled. We represent the solution 
of Eqs. 0 ,b) in the form 

X = Xie iult +X 2 e iU2t , <t> = 4>ie iLJlt + (/) 2 e iui2t , 0,b) 


where Xi 2 and <^ 12 are constant amplitudes. In the 
permanent oscillations regime, the reservoir population 
oscillates with a small amplitude around some average 
value N = (N) t . Further, we assume excitonic reservoir 


to be in the stationary state, N = 0. 

Let us first consider the phonon assisted pumping of 
the exciton component of the condensate. Substituting 
Eqs. (2a) and (|5j) into Eqs. & ,b) and separating real 
and complex parts we obtain 


RiN = Jx + 


ti 2 ! p 

7p/4 + ^1,2' 


w i,2 ~ ^ w i,2 — wi i2 (^ 2 - 7 p/4) - S^/p/A 


(6a) 
0. (6b) 


Clearly, only real roots of Eq. &) should be retained 
for the frequencies wi, 2 . 

Let us assume that the first requirement &) for the 
PT-symmetry is fulfilled, i.e., we set 6 = 0. In this case 
Eqs. (|6| admit equal amplitudes |A| = |xiI = |X 2 1 = 
10i | = | (^> 2 1 and characteristic polariton frequencies 


wi j2 = ±w = ± \Jfl 2 - 7p/4. (7) 

Equation ([ 7 ]) is familiar in PT-symmetry theory for cou¬ 
pled oscillators, H5HIZ1. The value of fl c = yp/2 can be 
associated with the PT-symmetry breaking threshold. In 
the case of a driven-dissipative exciton-polariton system, 
we operate significantly above the PT-threshold point, 
assuming that fl Q c , due to the strong exciton-photon 
coupling condition. Solution of Eq. (J6^i) with Eq. Q leads 
to the condition 


RiN = 7, (8) 

that simply implies the balance between pump and loss 
rates in the system and follows directly from the PT- 
symmetry criterion &)■ 

Figure [2] demonstrates the permanent oscillations 
regime for various pumping mechanisms of the excitonic 
component of the condensate. Permanent oscillations, 
shown in the left panel of the inset in Fig. [2] are estab¬ 
lished after several hundreds of picoseconds for exper¬ 
imentally accessible parameters. In particular, exciton 
and photon populations oscillate out of phase with equal 
constant amplitudes according to: |\| 2 = 4|A | 2 cos 2 (wf) 
and \cj >\ 2 = 4|A | 2 sin 2 (u;f), where |A | 2 = (P — P[ h ) / 2 7 
can be determined from Eqs. ([ljc) and | 8 |. Notably, the 
establishment of permanent Rabi oscillations regime is 
accompanied by the quick reservoir depletion (see the 
right panel on the inset in Fig. [ 2 ]) resulting in the sharp 
peak of photon/exciton populations for the first few pi¬ 
coseconds - see Fig. [2} Undoubtedly, since the conditions 
0 are satisfied now, this regime is also capable of hav¬ 
ing dynamical PT-symmetry properties occurring at the 
same time scale. 

For <5^0 Eq. (J^Jd) possesses three real roots, in gen¬ 
eral. Since the term S^p/A might be small enough for 
the moderate values of the detuning S we can omit this 
term and obtain (cf. 0 ) 

u,i >2 = 1 ± \J+ Ail 2 — 7 p^ . (9) 

The value of uip = \u>i — w 2 | represents the frequency of 
Rabi oscillations. However, since |aq| 7 ^ |w 2 |, Eq. (f(i|i) 
cannot be satisfied for both to i j2 simultaneously, cf. ul. 
Hence, we can conclude that permanent oscillations can¬ 
not be supported in the case of <5 7 ^ 0, which also corre¬ 
lates with the first PT- symmetry criterion (|4^). 

Now let us examine the permanent oscillations regime 
for the exciton pumping by exciton-exciton scattering. 
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FIG. 2. (Color online) The same as in Fig. [l]but for cw-pump 
P = 0.2 /xm -2 ps -1 , which is above threshold P th . The rele¬ 
vant behavior of reservoir is shown on the right panel of the 
inset. Permanent Rabi oscillations for excitonic and photonic 
components, plotted on a zoomed time scale, are shown on 
the left panel of the inset. 



In this case, using Eqs. © and omitting higher-order 
harmonics we obtain 

R2N 2 {\xi,2? + 2|x'2,l| 2 ) = lx + 2 !A Z P 2 ■ ( 10 ) 

7p/ z ± -h iu 12 

Notably Eq. (| 6 ]d) is valid in this case as well. Physically, 
it means that permanent Rabi oscillations frequency is 
independent on the physical mechanism of exciton pump¬ 
ing. The population N of the excitonic reservoir can be 
found from the equation 

P-'y R N-R 2 N 2 (|xt | 4 + |x 2 | 4 + 4|x 1 | 2 |x 2 | 2 ) = 0. (11) 


Again, in the limit of <5 = 0 we obtain a simple criterion 
3i? 2 A^ 2 |A | 2 = 7 of realization of permanent Rabi oscilla¬ 
tions, 

Figure [2] (yellow curve) demonstrates the population 
of excitonic \x\ 2 component in this case. The frequency 
cur of oscillations is still determined by Eq. 0. 

Now we consider permanent Rabi-oscillations in the 
case of non-zero detuning, { / 0 , and exciton pump¬ 
ing p 2 . Equations (10 1 admit solutions simultaneously 


for both u>\ and w 2 in the case of unequal amplitudes, 
Xi 7^ X' 2 - These solutions exist only for the absolute value 
of detuning S being smaller than some critical value 6 C . 
In fact, if d > |<5 C | population of an upper or a lower 
polariton branch (depending on the sign of 6) at the 
steady state becomes zero, that indicates the collapse 
of Rabi oscillations. It means that the initial system ([T]) 
has a steady state solution only; exciton y | 2 and photon 
\<j >\ 2 populations tend to some constant (non-zero) val¬ 
ues determined by this solution. The value of 5 C can be 
obtained using Eqs. ([9]) and (10) and it approximately 
reads 


fc=!± (iF 


lx - 8 ip + 3 ^/ 7 ^- + 81 p ^ 


( 12 ) 


FIG. 3. Time evolution of Stokes parameters for A z = 0.211, 
A = 0.1512. S x (red) and S y (blue) components - (a); S z 
(red) and So (blue) component - (b); (c) and (d) show the 
evolution of polarization on (S x , S^-plane (c) and on Stokes 
sphere (d) on the time scale of Tb = 2-k/ujb characterizing 
permanent Rabi oscillations beating period. The parameters 
are: px[N] = P 2 , P = 0.2 /rm _2 ps _ ; other values are the 
same as in Fig. [l] 


It is important that S c is determined by the parameters 
of a Rabi-oscillator as a whole, i.e., by the Rabi splitting 
frequency and damping rates jx, 7 p, and it does not 
depend on the parameters of the reservoir, including the 
external pump P. For GaAs-based semiconductor mi¬ 
crocavities we obtain HS C ~ 0.38/112 = 0.96 meV. Thus, 
permanent oscillations may be established in the vicinity 
of exciton-photon resonance for experimentally accessible 
parameters. 

Polarization properties. —Here we consider a com¬ 
plete set of Eqs. ([!]), that characterizes spin depen¬ 
dent exciton-polariton properties. Polarization proper¬ 
ties of the excitonic system are conveniently described 
by the Stokes vector S with components S XtVtZ defined 
by S = i (vErt . o-. ^r), where T = (x+,X-) T and a x>VjZ 
are the Pauli matrices. Since 7 ^ <5_ the behavior of 
polarization of exciton component is completely deter¬ 
mined by detuning A for a given non-zero Zeeman shift 
A z ¥=■ 0. 

The general case corresponding to both A and A z 
are nonzero is illustrated by Fig. [ 3 ] Using anzats © 
and definition for Stokes vector S it is easy to show that 
components S x and S y oscillate with the frequency A z 
(frequent beats in Figs. and [4^) . Besides there are 
fast oscillations with the frequency (ujr + + ojr-) /2 close 
to the frequency of Rabi-oscillations. If the conditions 
of establishing of permanent oscillations are fulfilled for 
both spin components simultaneously (|<5±| < S c ), the 
frequencies of Rabi oscillations lor+ and uir- are dif- 
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FIG. 4. Time evolution of Stokes parameters S x ^ ViZ for 
A z = A = 0.212. The parameters are: px[N] = Pi, 
P = 0.2 pm -2 ps -1 ; other values are the same as in Fig. Ill 


ferent, because of the relation |<5 + | ^ |J_|. In this case, 
the amplitudes of Stokes vector oscillations undergo addi¬ 
tional long-period beats with the period Tb = 2it/ujb = 
27r/ — ujr-) ~ see Fig.[3|b). These beats are 

sustained in time. 

In Fig. [3jd) we represented the evolution of the Stokes 
vector S normalised to the total number of excitons 


So = 


X+l 


lx-1 


on the Poincare sphere with a 

unit radius. Note that So, describing the excitonic sub¬ 
system, oscillates in time as shown in Fig. [3jb). 

Let us consider the limit where permanent oscilla¬ 
tions disappear for one of the spin components, i.e. for 
|<S + | < |i5 c | and |<5__| > |<5 C |, or vice versa, that can be 
achieved by the appropriate choice of A z and A. In 
fact, in this case the beats in Stokes parameters with the 
period Tb are suppressed in time and the amplitude of 
oscillations tends to the constant value. Figure ^demon¬ 
strates these features for Stokes parameters in the specific 
case of A = Az- In this limit, the detunings 5- = 0 and 
(5+ = 2A^. Physically, it means that permanent oscilla¬ 
tions can be established for an exciton pumping mecha¬ 
nism described by (2a), but only for ” spin component 
while population of another component will approach a 
constant value. At long time intervals only fast oscilla¬ 
tions survive. 

In conclusion, we have described the regime of per¬ 
manent Rabi oscillations in a driven-dissipative exciton- 
photon system feeded from the incoherent exciton reser¬ 
voir. Two types of reservoir-Rabi-oscillator coupling 
mechanisms have been examined. We have shown that 
the permanent Rabi oscillation regime is established 
above the threshold pumping where the dynamical PT- 
symmetry of the coupled exciton-photon system is re¬ 
alised. Notably, this regime occurs in a non-resonant 
case only if the exciton component of the condensate 
is pumped by exciton-exciton scattering predominantly. 
We studied excitonic polarization properties in the pres¬ 
ence of external magnetic fields. The conditions for per¬ 
manent Rabi oscillations for polarized excitonic and/or 
photonic system are established and discussed. Realiza¬ 
tion of permanent Rabi oscillations in this case play a cru¬ 
cial role for creation and manipulating of long-lived spin 
polarisation in the exciton-photon system which pave the 
way to the design of new optical memory devices. 
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